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ABSTRACT 


A  derivation  of  the  characteristic  equations  for  supersonic  inviscid 
flow  with  magnetohydrodynamic  (MHD)  forces  present  is  given  for  the 
two-dimensional  and  axisymmetric  cases.  Workable  forms  of  the  equa¬ 
tions  relating  the  MHD  effects  are  indicated  for  the  axisymmetric 
problem. 
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a  Sonic  speed 

B  Angular  component  of  B 

B  Magnetic  induction 

E  Electric  intensity 

Ex  x  component  of  E 

— ► 

Ey  y  component  of  E 

H  Total  enthalpy 

h  Enthalpy 

Current  density 

— > 

Jx  x  component  of  J 

— ► 

Jy  y  component  of  J 
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p  Pressure 

s  Entropy 

u  x  component  of  v 

v  y  component  of  v 

v  Velocity 

x  Axial  coordinate 

y  Radial  coordinate 

Pq  Permeability  of  free  space 

p  Density 

a  Conductivity 

4>  Set  of  variables  (x,  y,  u,  v,  p,  and  h) 
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SECTION  I 
INTRODUCTION 


The  conservation  equations  for  the  axisymmetric  supersonic  inviscid 
flow  problem  with  magnetohydrodynamic  (MHD)  forces  present  are  taken 
to  be 


dp  \ 

c  v  T  v  — —  + 
•  3y 


3  u  3  u 

P  u  ; —  +  P  v  - — 

d  x  ay 


d  \  d  \ 

p  u  -  +  p  v  - 

dx  dy 


+ 


t£JL  =  0 

$  X 


(1) 

(2) 

(3) 


,u|iL  +  =  j 

<7  x  d  y 


E 


(4) 


The  primary  purpose  of  the  work  that  follows  is  deriving  the  char¬ 
acteristic  equations  corresponding  to  the  above  system  and  noting  an 
appropriate  set  of  relations  for  the  MHD  effects  that  would  be  solved 
simultaneously  with  the  characteristic  equations. 

The  manner  in  which  the  results  of  the  following  derivation  can  be 
made  applicable  to  the  two-dimensional  problem  is  indicated  in 
Section  V. 


Coupled  with  the  conservation  equations  are  five  auxiliary  relations 
which  must  also  hold. 


h  +  T 

■  (u2  +  V2) 

(5) 

P  =  fo 

(  p,  h  ) 

(6) 

J>B  - 

fi  (f) 

(7) 

JxB  = 

*2  V<£) 

(8) 

— »  -4 

J  .  E  = 

U  (4>) 

(9) 

Equation  (5)  is  the  definition  of  total  enthalpy.  Equation  (6)  is  the 
equation  of  state.  It  is  assumed  that  fg  has  continuous  first  derivatives 
with  respect  to  p  and  h.  Equations  (7),  (8),  and  (9)  are  intended  to 
assert  that  the  MHD  forces  are  known  continuous  functions  of  d. 
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The  reader  is  reminded  that,  for  the  axisymmetric  problem,  f]_,  f2, 
and  f3  must  obey  the  appropriate  laws  of  symmetry.  1 

An  explicit  set  of  conditions  for  the  MHD  relations  is  noted  in 
Section  VI.  For  the  derivation  of  the  characteristic  equations,  however, 
it  is  convenient  not  to  introduce  these  cumbersome  expressions. 

As  this  is  a  problem  of  deriving  the  characteristic  equations  cor¬ 
responding  to  Eqs.  (1),  (2),  (3),  and  (4)  and  not  one  of  actually  solving 
the  equations,  boundary  and  initial  conditions  need  not  be  specified. 


SECTION  II 

CHOICE  OF  VARIABLES 


The  independent  variables  in  Eqs.  (1),  (2),  (3),  and  (4)  are  x  and  y. 
The  dependent  variables  are  u,  v,  p,  and  h,  partly  as  a  matter  of  choice. 
It  is  desirable,  therefore,  to  rewrite  the  system  of  partial  differential 
equations  so  that  derivatives  with  respect  to  the  independent  variables  of 
only  the  dependent  variables  appear.  For  example,  since  p  is  not  a 
dependent  variable  and  in  lieu  of  Eq.  (6) 

dp  dp  dp  d  p  d h 
d  x  dp  d  x  d  h  d  x 


and  can  he  obtained  from  Eq,  (6),  and  and  are  derivatives 
dp  dh  dx  9x 

of  dependent  variables  with  respect  to  an  independent  variable  (dependent 

derivatives). 


The  desired  form  of  Eqs.  (1),  (2),  (3),  and  (4)  is‘ 


d u  d v  d  o  d p  d  p  d p  d  p  d  h 

py  -  +  py  -  +  uy  - —  +  vv  +  yu  — — 


dy 


dh  d: 


dp  d  h 

+  vv  — —  - —  =  —  pv 
d  h  d  v 


d u  d  u  dp  ID 

P  u  - —  +  pv  —  +  — r-  =  J  B 

a  x  dy  d  x  1 

d \  d  v  dp  , 

pu  +  pv  +  — =  -JXB 

ax  d y  d  y 


du 


d  v 


dh 


dh 


P —  +  puv  ^  +  puv  +  pvJ  +  pu  +  pv  2JL  =  j  .  E 

d  x  dy  dx  dy  dx  dy 


(10) 

(11) 

(12) 

(13) 


^i.  e.  ,  fi  and  f3  are  even  functions  of  y,  and  f2  is  odd. 

^It  is  assumed  until  otherwise  indicated  that  y  ^  0.  A  limiting 
process  will  he  employed  in  Section  IV,  Part  III  to  determine  certain 
conditions  on  the  x-axis. 
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Notice  that  the  above  equations  can  be  interpreted  as  a  system  of 
four  linear  algebraic  equations  in  eight  unknowns,  the  unknowns  being 
the  dependent  derivatives.  Also,  by  definition. 


dx  T 
9  k 

i-2-  dv  =  du 
d  y 

(14) 

i-2-  dv  + 
d  \ 

dv  =  d  V 

d  > 

(15) 

dP  , 

— —  d  v  + 
d  x 

4^  dv  =  dp 

dy 

(16) 

ii!-  d*  - 

dv  =  dh 

(17) 

d  x 

dy 

whenever  these  differentials  exist.  Equations  (10)  through  (17)  are, 
therefore,  a  system  of  eight  linear  algebraic  equations  in  eight  un¬ 
knowns  . 

This  concept,  that  of  coupling  the  given  differential  equations  with 
-the  equations  defining  the  dependent  differentials,  and  interpreting  the 
results  as  a  system  of  linear  algebraic  equations,  is  fundamental  to 
this  method  of  deriving  the  characteristic  equations. 

As  much  of  what  follows  is  directly  related  to  the  system  of  eight 
equations,  they  are  presented  according  to  a  convenient  format  in  the 
following  illustration.  ^ 


f^U 

fix' 

fill 

ay 

<K 

C’  V 

0p  |  dp 

?)x  J  r>v 

c)h 

fix 

■Bh 

ay 

py 

0 

_ 

n 

py 

rip 

5p 

J  9p 

y  u  — 
r*li 

>V9h“ 

-pv 

pU 

.  Pv  o 

0  '  ] 

0 

0 

0 

JjB 

0 

0 

pu  ]  pv  I  0 

1  1  0 

t 

0 

-,rxB 

9 

p  u  p  uv 

pu-.-  |  pv2 

0  1  0  pu 

pv 

i  .  u 

dx 

dy 

0 

0 

0 

0 

0 

0 

du 

0 

0 

<\K 

dv 

0 

0 

0 

0  . 

d\ 

0 

0 

0 

0 

dx 

d.v 

0 

0 

dp 

0  0  0  0  0 

L  1  .  _ 1 

0 

dx  d\ 

■  ‘ 

dh 

3 


Shapiro  made  use  of  such  a  format  in  Ref, 


1,  p.  518. 
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SECTION  III 

DEFINITION  OF  THE  CHARACTERISTICS 


With  the  introduction  of  a  limited  amount  of  matrix  notation,  the 
definition  of  the  characteristic  equations  can  be  conveniently  and  simply 
stated. 

— ^ 

Let  X  be  the  8  by  1  colurmqvector,  the  components  of  which  are  the 
dependent  derivatives,  and  let  Y  be  the  8  by  1  column  vector  whose 
elements  are  the  right-hand  side  of  the  system  of  equations, 


Let  A  be  the  8  by  8  coefficient  matrix.  Then  Eqs.  (10)  through  (17)  can 
simply  be  written 

A  X  =  Y  (18) 

The  A  augmented  matrix  will  be  denoted  by  C4  That  is,  C  is  the  8  by  9 
matrix  obtained  by  including  the  elements  of  Y  as  an  additional  column 
to  the  matrix  A. 

Definition1^:  The  set  of  characteristic  equations  consists  of  all 
equations  that  follow  on  requiring  that  R(A)  =  R(C)  <  7. 

This  definition  is  motivated  by  the  question,  "Do  there  exist  curves 
within  the  domain  of  interest  on  which  Eq.  (18)  has  multiple  solutions  ?  " 
It  is  known  from  the  theory  of  linear  algebra  that  if  such  curves  exist, 
then  it  is  necessary  that  R(A)  £  7  (Ref.  2,  p.  61).  It  is  known  also,  that 
no  solution  can  exist  unless  R(A)  =  R(C)  (Ref.  3,  p.  15). 


SECTION  IV 

DERIVATION  OF  THE  CHARACTERISTIC  EQUATIONS 


The  purpose  of  this  section  is  to  determine  the  equations  that  follow 
on  requiring  that  R(A)  =  R(C)  <  7.  These  results  will  be  derived  in 


^For  any  matrix  K,  R(K)  denotes  the  rank  of  K. 
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Part  I  and  Part  II.  Certain  special  cases  will  be  treated  in  Part  III. 


4.1  PART  I 


Since  C  is  the  A  augmented  matrix,  it  follows  that  R(A)  <  R(C) 
necessarily  (Ref.  3,  p.  15). 

Let  B  be  the  8  by  8  matrix  obtained  by  replacing  any  column  of  A 
(say  the  first  column  for  definiteness)  with  the  elements  of  Y.  Then 
R(C)  £  7  only  if  det6  (A)  -  0  and  det  (B)  =  0  simultaneously. 

The  result  of  setting  det  (A)  =  0  and  simplifying  is 

(udy  -  vdx)2  [( u2  -  a*)dy2  -  2uvdydx  -  <  v*  -  a1)  dx*]  =  0 

where,  because  of  the  choice  of  dependent  variables6, 

a2  =  - P _ 

dP  dp 

<H  d  P 

From  Eq.  (19)  either’ 

□dy  -  vdx  -  0  (21) 

or 

(u2  -  a’)  dy*  -  2uvdyd\  +  (v2  -  a2)  d\2  -  0  (22) 


(19) 


(20) 


Equation  (21)  is  the  equation  of  a  streamline.  Equation  (22),  a 
quadratic  in  dy/dx,  is  commonly  referred  to  as  the  equation  of  a  left  or 
right  running  characteristic.  Setting  det  (B)  -  0  gives 


f  udy  —  vdx)  U  j7j^Bdv  —  pvdu )  (udv  — 


^dydpj 


+  4*  dy  (JxBvdy  +  J  ■  E  dy  -  J  Budv)  +  -^-1  dv2( 

dh  >  y  \ 


f  pdy(dxdp  +  pudydv  +  JxBdydx) 
-  pdx  fJyBdv  -  pvdu)  =  0 


(23) 


6The  determinant  of  A 

6Appendix  I  gives  a  proof  that  (a)  as  defined  by  Eq.  (20)  is  the  speed 
of  sound  provided  additional  assumptions  on  the  nature  of  the  fluid  are 
made. 
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Since  det  (A)  =  0  and  det  (B)  =  0  must  hold  simultaneously,  the  form  of 
Eq.  (23)  can  be  simplified  somewhat  by  imposing  Eq.  (22)  as  a  side 
condition. 

( udy  -  \dx)  ( udv  -  vdu)  +  ( udy  -  vdx)2  — 

r. , 

4.  -iP  (vdv  +  udx)  +  -U  (vv@+  j  ■  E  -  uj  B)  (udy  -  vdx)2 

p  '  •  p  dy  ’x  > 

+  -5-  ( J  dx  +  J  dv)  (udv  -  vdx)  =  0  (24) 

P  \  »  ' 

Equation  (24)  is  known  as  the  compatibility  equation  along  a  left  or  right 
running  characteristic. 

When  Eqs.  (22)  and  (24)  are  satisfied  simultaneously,  R(B)  =  R(A)  -  7 
by  inspection  and  R(C)  =  7  (Ref.  4,  p.  54). 


4.2  PART  II 


~\ 


Along  a  streamline  (when  udy  -vdx  =  0)  R(A)  =  6  by  inspection,  so 
that  R(C)  is  at  most  seven.  In  this  case  it  can  be  shown  that  R(C)  =  6  if 
and  only  if 

-if  -  udu  +  vdv  +  (JxBdy  -  JyBdx) j  p  =  0  (25) 


and 


dh  +  udu  +  vdv  =  J-V  f  Ai. 


(26) 


Satisfying  Eqs.  (25)  and  (26)  does  not  affect  R(A);  hence  Eqs.  (24),  (25), 
and  (26)  are  the  desired  relations  on  the  streamline. 


4.3  PART  III 

The  characteristic  equations  for  the  two-dimensional  case  are 
exactly  the  same  as  those  previously  derived  except  that  the  term  con¬ 
taining  the  factor  v/y  in  Eq.  (24)  is  replaced  with  zero. 

For  the  axisymmetric  case,  Eq.  (24)  is  singular  on  the  x  axis 
because  of  the  v/y  term.  Applying  THospitaFs  rule  and  making  use  of 
Eq.  (1),  it  can  be  shown  that 


/' 

/ 
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lim  — 


I  dp  ii. 


y  -*  0  y  2  pd 


Since  p  is  not  a  dependent  variable,  a  more  workable  form  is 

lim  -  =  -i-  V—  -  —  -  -  /uj  B  ^  +  LI  ifM 

v-°  >'  2  i\a*  I  dx  P  \  ■  *  P  P  ^h/(y  =  1 


Equation  (28)  was  obtained  by  differentiating  the  product  in  Eq.  (27)  and 
utilizing  Eqs.  (6),  (20),  (25),  and  (26)  (note  that  y  =  0  is  a  streamline). 
Define 


-f  if  v  *  0 


lim  —  if  y  =  0  ( 
v^o  y  } 


(  axisy 


symmetric  case 


}  two-dimensional  case 


The  characteristic  equations  for  either  the  axisymmetric  or  two- 
dimensional  case  are  now  summarized  for  convenience.  Along  a  left  or 
right  running  characteristic 


(u2  -  a2)  dvJ  -  2uvdydx  -  (  y2  —  a2)  dx1  =  0 


(udy  -  vdx)(udv  -  vdu)  +  lady  -  vdx)2fC 


+  — -  (vdy  +  udx )  +  —5—  (vj  B  +  J  •  E  -  uLB)  ( udv  -  vdx)2.  c/\ 

p  p2  d  h  31  y  -  \  '  J 


iT  au  rL 


+  —  (Jxdx  +  J  ydy )  ( u  dy  -  vdx)  =  0 


and  along  a  streamline 


udy  -  vdx  -  0  SjanJi'i  dfr§tWux->@nAJt  (32) 

L  .  udu  +  vdv  +  UxBdy  -  JvBdx)  /  p  =  0  (33) 


<!  IWify  -Cj  tot  h’&y  •  dh  +  udu 


vdv  = 


J  •  E  d  x 
P  u 


An  independent  check  of  the  above  compatibility  equations  is  given 
in  Appendix  II. 
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SECTION  V 

MHO  RELATIONS  FOR  THE  AXISYMMETRIC  PROBLEM 


Maxwell's  equation  along  with  Ohm's  law  are  the  only  relations 
needed  to  obtain  the  MHD  terms  appearing  in  Eqs.  (31),  (33),  and  (34) 
assuming  adequate  boundary  conditions  are  given  and  assuming  that 
a  is  a  known  function  of 


V.  E 

=  0 

(35) 

V  x  E 

=  0 

(36) 

V-  B 

=  0 

(37) 

VxB 

=  he  J 

(38) 

J  =  0 

(E  +  V  x  B) 

(39) 

Equations  (35)  and  (36),  an  elliptic  system  of  equations,  can  be  un¬ 
coupled  from  the  remaining  equations.  If  "reasonable"  boundary  con¬ 
ditions  are  imposed  on  these  two  equations,  a  closed  form  solution  can 
be  obtained  giving  Ex  and  Ey  as  functions  of  position  only.  Hence, 
assuming  E  to  be  a  known  function  of  position,  an  attempt  can  be  made 
to  solve  the  scalar  equations  corresponding  to  Eqs.  (37),  (38),  and  (39) 
simultaneously  with  the  system  of  characteristic  equations  (Eqs,  (30) 
through  (34)). 


Equation  (37)  holds  identically,  and  Eqs.  (38)  and  (39)  can  be  written 


d  y  B  , 

ST  =  -Wty 

(40) 

d  y  B 

By  = 

(41) 

J,  =  <r(Ex  +  vB) 

(42) 

Jy  =  a(Ky  -  uB) 

(43) 

of  MHD  relations  can  be  summarized  as 

°  =  gj  (<£) 

(44) 

Ex  =  g2  U,  >') 

(45) 

Ey  =  (  x>  y  ) 

(46) 

Jx  =  £X(EX  +  vB) 

(47) 

Jy  =  O  ( Ev  -  uB) 

(48) 

=  |(Ex  +  vB)  ydy  -  (Ey  - 

uB)  ydx  | 

(49) 
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Equation  (49)  was  obtained  by  using  Eqs.  (40)  and  (41)  to  obtain  the  total 
differential  of  yB  and  by  using  Eqs.  (42)  and  (43)  to  eliminate  Jx  and  Jy. 
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APPENDIX  I 

SONIC  SPEED  RELATIONS 


It  was  mentioned  earlier  that  the  only  assumption  necessary  on  the 
type  of  gas  being  considered  is  that  p  be  a  given  function  of  p  and  h  with 
continuous  first  derivatives. 


P  =  M  p,  li ) 

In  order  to  prove  that  (a)  as  given  by  the  relation 


*£-\  +  Pi£.\ 


(i-i) 


(1-2) 


is  the  sonic  speed,  certain  additional  assumptions  are  made. 


It  is  assumed  that  a  differentiable  relation  defining  entropy  exists 

s  =  ea  (P-  h)  (1-3) 

and  that  unique  inverse  relations  giving  p  and  h  as  functions  of  p  and  s 
also  exist.  The  definitions  of  p  and  a  are  taken  to  be 

o-lf)  H-« 

»'  -  £r),  (I‘5> 

Finally,  it  is  assumed  that 


It  will  now'  be  shown  that  the  terms  on  the  right-hand  side  of 
Eqs.  (1-2)  and  (1-5)  are  equivalent. 


Differentiating  Eqs.  (1-1)  and  (1-3)  with  respect  to  p  at  constant  s 
yields 


=  0 


dk\  _dp\ 

Ms  Ms 


=  1 
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Since  the  Jacobian  related  to  the  above  equations  cannot  vanish,  it 
follows  that 


This  completes  the  proof. 
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APPENDIX  II 

CHECK  ON  THE  COMPATIBILITY  EQUATIONS 


Because  of  the  involved  process  of  obtaining  the  characteristic  equa¬ 
tions,  an  independent  and  relatively  simple  check  of  the  compatibility 
equations  is  offered. 

Consider  Eq.  (34)  which,  according  to  a  previous  assertion,  is  a 
compatibility  equation  on  a  streamline.  This  equation  can  be  checked  as 
follows : 

1.  Multiply  Eq.  (4)  through  by  dx. 

2.  Impose  the  condition  udy  =  vdx  (Eq.  (32))  on  the  result. 

3.  Note  the  definition  of  total  enthalpy  (Eq.  (5)). 

4.  Equation  (34)  follows. 

Equation  (33),  the  other  compatibility  equation  on  a  streamline,  can 
be  checked  in  much  the  same  way: 

1.  Multiply  Eqs.  (2)  and  (3)  by  dx  and  dy,  respectively. 

2.  Impose  the  condition  udy  -  vdx  on  each  equation. 

3.  Add  the  resulting  equations. 

4.  Equation  (33)  follows. 


The  check  on  the  compatibility  equation  along  a  left  or  right  runnin 
characteristic  (Eq.  (31))  is  somewhat  more  involved: 


1. 

Multiply  Eq.  (1)  by  (udy  -  vdx)2,  (  1  )  . 

2. 

Multiply  Eq.  (2)  by  j 

—  (udy  -  vdx)2  -  dy  (udy  -  vdx) 

p  on 

1 

TT 

3. 

Multiply  Eq.  (3)  by 

r  g 

|v/p  (udy  -  vdx)2  +  dx  (udy  -  vdx) 

4. 

Multiply  Eq.  (4)  by 

[■  p  9h  (udy  '  Vdx)2]  P  ’ 

g 


v  .  Add  the  results  of  items  1  through  4  above  and  simplify. 

6.  Impose  Eq.  (30)  as  a  side  condition  on  the  result  of  the  above. 

7.  Equation  (31)  follows. 
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